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We aim at understanding transport in porous materials including regions with both high 
and low difhisivities. For such scenarios, the transport becomes structured (here: micro- 
macro). The geometry we have in mind includes regions of low diffusivity arranged in a 
locally-periodic fashion. We choose a prototypical advection-diffusion system (of minimal 
size), discuss its formal homogenization (the heterogenous medium being now assumed to 
be made of zones with circular areas of low diffusivity of x- varying sizes), and prove the 
weak solvability of the limit two-scale reaction-diffusion model. A special feature of our 
analysis is that most of the basic estimates (positivity, L°°-bounds, uniqueness, energy 
inequality) are obtained in ^-dependent Bochner spaces. 
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1 Introduction 

We consider transport in heterogeneous media presenting regions with high and low 
diffusivities. Examples of such media are concrete and scavenger packaging materials. 
For the scenario we have in mind, the old classical idea to replace the heterogeneous 
medium by a homogeneous equivalent representation (see [TJ [2j [5j [22] and references 
therein) that gives the average behaviour of the medium submitted to a macroscopic 
boundary condition is not working anymore. Specifically, now the transport becomes 
structured (here: micro-macr^ [3] |14| . 

The geometry we have in mind includes space-dependent perforation^] arranged in a 
locally-periodic fashion. We refer the reader to section [2] (in particular to Fig. [T]) , where 
we explain our concept of local periodicity. Our approach is based on the one developed 
in 24 25 and is conceptually related to, e.g., [6 11 . When periodicity is lacking, the 



typical strategy would be to tackle the matter from the percolation theory perspective 

1 "Micro" refers here to a continuum description of a porous subdomain at a separated (lower) 
spatial scale compared to the "macro" one. 

2 By "space-dependent perforations", we mean that at each spatial position x, our model will 
allow us to zoom in a a;-dependent pore space, or subject to a more general interpretation, a 
^-dependent porous subdomain, called here perforation. 
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(see e.g. chapter 2 in 12 and references cited therein^]) or to reformulate the oscillating 
problem in terms of stochastic homogenization (see e.g. [I]). In this paper, we stay within 
a deterministic framework by deviating in a controlled manner (made precise in section 
[2| from the purely periodic homogenization. 

We show our working methodology for a prototypical diffusion system of minimal size. 
To keep presentation simple, our scenario does not include chemistry. With minimal 
effort, both our asymptotic technique and analysis can be extended to account for vol- 
ume and surface reaction production terms and other linear micro-macro transmission 
conditions. We only emphasize the fact that if chemical reactions take place, then most 
likely that they will be hosted by the micro-structures of the low-diffusivity regions. We 
discuss the microscale model for the particular case in which the heterogenous medium 
is only composed of zones with circular areas of low diffusivity of x-varying sizes. This 
assumption on the geometry should not be seen as a restriction. We only use it for ease 
of presentations and it does not play a role in our formal and analytical results. Our 
asymptotic strategy is based on a suitable expansion (remotely resembling the boundary 
unfolding operator [7]) of the boundary of the perforations in terms of level-set func- 
tions. In particular, we can treat in a quite similar way situations when free-interfaces 
travel the microstructure; we refer the reader to [24| for a dissolution precipitation free- 



boundary problem and 20 for a fast-reaction slow-diffusion scenario where we addressed 
the matter. 

The results or our paper are twofold: 

(i) We develop a strategy to deal (formally) with the asymptotics e — > for a locally 
periodic medium (where e > is the microstructure width) and derive a macroscopic 
equation and x-dependent effective transport coefficients (porosity, permeability, tor- 
tuosity) for the species undergoing fast transport (i.e. that one living in high diffusivity 
areas), while we preserve the precise geometry of the microstructure and correspond- 
ing balance equation. The result of this homogenization procedure is a distributcd- 
microstructure model in the terminology of R. E. Showalter, which we refer here as 
two-scale model. 

(ii) We analyze the solvability of the resulting two-scale model. Solutions of the two- 
scale model are elements of re-dependent Bochner spaces. Our approach benefits from 
previous work on two-scale models by, e.g., Showalter and Walkington [23], Eck [9], 
and Meier and Bohm [17 18 . A special feature of our analysis is that most of the 



basic estimates (positivity, L°°-bounds, uniqueness, energy inequality) are obtained in 
the x-dependent Bochner spaces. Our existence proof is constructed using a Schaudcr 



fixed-point argument and is an alternative to 23 , where the situation is formulated as 
a Cauchy problem in Hilbert spaces and then resolved by holomorphic semigroups, or 
to 17 , where a Banach-fixed point argument for the problem stated in transformed 



domains (i.e. x-independent) is employed. 
Note that (i) and (ii) are preliminary results preparing the framework for rigorously 



3 Fig. 2.3 (a) from 12 , p. 39 illustrates a computer simulation of the consolidation of spherical 
grains showing regions with high and low porosities corresponding to high and low diffusivity 
areas. 
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Figure 1. Schematic representation of a locally-periodic heterogeneous medium. The cen- 
ters of the gray circles are on a grid with width e. These circles represent the areas of 
low diffusivity and their radii may vary. 

proving a convergence rate for the asymptotics e — > 0; we will address this convergence 
issue elsewhere. 

The paper is organized in the following fashion: Section [2] contains the description 
of the model equations at the micro scale together with the precise geometry of our 
x-dependent microstructure for the particular case of circular perforations. The homog- 
enization procedure is detailed in section [3] The main result of this part of the paper 
is the two-scale model equations as well as a couple of effective coefficients reported in 
section [4j The second part of the paper focusses on the analysis of the two-scale model; 
see section [5] The main result, i.e. Theorem 5.11 ensures the global-in-time existence 



of weak solutions to our two-scale model and appears at the end of section 5.3 A brief 
discussion section concludes the paper. 



2 Model equations 

We consider a heterogenous medium consisting of areas of high and low diffusivity. The 
medium is in the present paper represented by a two dimensional domain. We denote 
the two dimensional bounded domain by fi C K 2 , with boundary T, and for ease of 
presentation we suppose in this section that the areas of the medium with low diffusivity 
are circles. We do not use this restriction in later sections; there the areas with low 
diffusivity can have different shapes, as long as neighboring areas do not touch each 
other. 

Let the centers of the circles with low diffusivity, with radius Rij < e/2, be located 
in a equidistant grid with nodes at (ei, ej), where e is a small dimensionless length scale. 
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We assume that there is given a function r(x) : O — >• [0, 1/2) such that the radii Rij of 
the circles Bij are given by Rij = er(xij), where Xij = (ei,ej). We define the area of 
low diffusivity Q e , which is the collection of the circles of low diffusivity, as £l\ := UBij 
and we define the area of high diffusivity fi^, which is the complement of flf in f2, as 
Vi e h := Q\f2 e . The boundary between high and low diffusivity areas is denoted by T e , 
which is given by T f := d£lf. It is important to note that we assume that the circles of 
low diffusivity do not touch each other, so that n Tki = if i ^ k or j ^ I, where 
Tij := dBij, and we also assume that the area of low permeability does not intersect the 
outer boundary of the domain fi, so that T n = for all 

We denote the tracer concentration in the high diffusivity area by u € , the concentration 
in the low diffusivity area by v e , the velocity of the fluid phase by q e and the pressure by 
p e . All these unknowns are dimensionless. In the high diffusivity area we assume for the 
fluid flow a Darcy-like law and incompressibility, while we neglect fluid flow in the low 
diffusivity area. The diffusion coefficient in the low diffusivity area is assumed to be of 
the order of 0(e 2 ), while all the remaining coefficients are of the order of 0(1) in e. We 
assume continuity of concentration and of fluxes across the boundary between the high 
and low diffusivity areas. 

The model is now given by 



V • (D h Vu £ - q e u e ) 

(2.1) 
(2.2) 
(2.3) 

(2.4) 
(2.5) 



q e — —nVp e 
V • q £ = 


in 


v\ = e 2 V • (D,Vd £ ) 


in ttf, 


v f - ■ (D h Wu e ) = eV • 
w £ = v e 


(AV« £ ) 

on r £ 


q e = 




u e (x, t) — Ub(x, t) 
q e (x,t) = q b (x,t) 


on r, 


u e (x, 0) = u e j(x) in 
v e (x, 0) = vj(x) in 





where denotes the diffusion coefficient in the high diffusivity region, D\ the diffusion 
coefficient in the low diffusivity regions, k denotes the permeability in the Darcy law for 
the flow in the high diffusivity region, v e denotes the unit normal to the boundary T e (t) , 
where q b and U], denote the Dirichlct boundary data for the concentration u e and Darcy 
velocity q e and where u e j and vj denote initial value data for the concentration u e and 
v e . 
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3 Formal homogenization 

In addition to the macroscopic variable x, we introduce a periodic unit cube U with 
microscopic variable y: 

y = (y u y 2 ), and U := {y e K 2 | - 1/2 < Vi < 1/2 for % = 1,2}). (3.1) 

For the formal homogenization we assume the following formal asymptotic expansions 
for u e , v e , q e and p t : 

u £ (x, t) — uo(x. x/e, t) + eu\(x, x/e, t) + e 2 u 2 (x, x/e, t) + ... 
v e (x, t) — V(s{x, x/e, t) + evi(x, x/e, t) + e 2 vi{x, x/e, t) + ... 
q e (x,t) = q (x,x/e,t) + eqi(x,x/e,t) + e 2 q 2 {x,x/e,t) + ... 
p e (x,t) = po(x,x/e,t) + epi{x,x/e,t) + e 2 p 2 {x, x/e, t) + ... 

where u k {-,y, ■), v k (-,y, ■), q k {-,y, •) and Pk(-,y, •) are 1-periodic in y = f . The gradient 
of a function f(x, depending on x and y = ~ is given by 

V/ = V a! / + iv»/| w= «, (3.2) 

where and V y denote the gradients with respect to the first and second variables of 
/■ 



3.1 Level set formulation of the perforations boundary 

Since the location of the interfaces between the low and the high diffusivity regions also 
depends on e, we need an e-dependent parametrization of these interfaces. A convenient 
way to parameterize the interfaces is to use a level set function, which we denote by 
S*(x): ' 

ier e » s e (x) = o. 

Since we allow the size and shape of the perforations to vary with the macroscopic variable 
x, we might use the following characterization of S e : 

S e (x) = S(x,x/e) (3.3) 

where S : O X U — > K is 1-periodic in its second variable, and is independent of e. In 
this section we show, using the example of a grid of circles with varying sizes, that this 
characterization of S* e is not sufficient to characterize all locally-periodic sequences of 
perforation geometries. In fact, we need to expand S e as 

S e (x) = S Q (x,x/e) +eS 1 {x,x/e) + e 2 S 2 (x,x/e) + ... (3.4) 

where the Si : 57 x U —> K are 1-periodic in their second variable, for i = 0,1, 2, ... and 
are independent of e. 

In order to find an explicit expression for S e (x) in this particular case, i.e. the case of 
circular domains with radius r(x) (see Fig.[T]), we define P(x) to be the periodic extension 
of the function x — > \x\ and Q{x) to be the periodic extension of the function x — > x, 
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both defined on the square [— i , A] x [— A, |], given by 

P{x) = P( Xl ,x 2 ) = VL^i + 1/2J 2 + L^2 + 1/2J 2 , 
Q(x) = Q(a;i,a; 2 ) = ([an + 1/2J, [a* + 1/2J), 

where [aj := max{n G Z | n < a} denotes the floor of a (rounding down). We can write 
S^x) as follows: 

S e (a;) =r(x-eQ(x/e))-P(x/e). (3.5) 



Interestingly, the expression (3.5) plays the same role as the boundary unfolding operator 
(cf., for instance, [7| Definition 5.1). Note that S e is not a continuous function, it jumps 
when x\ or x 2 cross a multiple of e. Whenever we assume that r(x,t) < 1/2, this is not 
a problem, since in this case S e is continuous and smooth in a neighborhood of its zero 
level set, which is what we are interested in. 

To check that the zero level set of S e consists indeed of circles around X{j with radius 
er(xij), we consider a curve, which without loss of generality can be parametrized in the 
square with sides e around Xij by Xij +7(3). For this curve to be a zero level set, it should 
hold that 

r(xij + 7(8) - eQ{e-\ XlJ + 7 (s)))) = Pfc-^sy + 7 (s))). 
Using that x^ — (ei, ej), with ei, ej G eZ (~l f2, we obtain 

r((a, ei) + 7(a) - cQ((*,i) + e _1 7(s))) = P{{h j) + 7(*)). 
and using the periodicity of P and Q we get 

= ItO)], 

which means that 7 (s) should be a circle with radius r(xjj). 

Now we can write the level set function S e formally as the expansion 

S e {x) = S {x,x/e) + eSi(x,x/e) + e 2 S 2 (x,x/e) + (9(e 3 ), 

where Sk(-, y, for k = 0, 1, 2, are 1-periodic in y = and are independent of e. In 
order to find the terms in this expansion, we assume that r is sufficiently smooth and so 
that we can use the Taylor series of r around x: 

r(x - eQ(x/e) = r(x) - eQ(x/e) ■ Vr(i) + —Q(x/e) ■ V 2 r(x)Q(x/e) + 0(e 3 ), 

where V 2 r denotes the Hessian of r w.r.t. x. This suggests the following definition of the 
terms in the expansion of S e : 

S (x,x/e) := r(x) - P{x/e), 
Si(x, x/e) := —Q(x/e) ■ Vr(x), 

S 2 (x,x/e) := ~Q(x/e)-V 2 r(x)Q(x/e), 
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In ( 2.3 l ) we have used the superscript e for the normal vector v e in the interface conditions 
for v e and u e . The reason is that the normal vector depends on the geometry of the 
different regions, and this in turn depends on e. In order to perform the steps of formal 
homogenization, we have to expand v e in a power series in e. This can be done in terms 
of the level set function S e : 

S= Iiy' X/ ,i atxel". (3-6) 
\VS'{x iX /e)\ 



First we expand |VS e |. Using the chain rule (3.2) (see also [12]), the expansion of S e and 
the Taylor series of the square-root function, we obtain 

\VS^\ = - e \V y So\ + 0(e°). (3.7) 

In the same fashion, we get 

v € = v Q + evi + 0(e 2 ), 

and 

VzSq + Vj,Si (V^^o • VyS'o + VyS'o • Vj,5i) V^So 



where 



|V y S„| |V,5 P |V,5o|' 

If we introduce the normalized tangential vector To, with To _L v$, we can rewrite v\ as 

Vl ~ T ° IV^oI ' (3 ' 8) 

Now we focus on the interface conditions posed at T e . In order to obtain interface condi- 
tions in the auxiliary problems, we substitute the expansions of u e , q t , and u e into (2.3). 
This is not so straight-forward as it may seem, since the interface conditions ( |2.3 1 arc 
enforced at the oscillating interface T e , i.e. at every x where S e (x) = 0. For formulating 
the upscaled model it would be convenient to have boundary conditions enforced at 

r Q (x):={y\S (x,y)=0}. (3.9) 

To obtain them, we suppose that we can parametrize the part of the boundary r|. that 
surrounds the sphere Bij with fc £ (s), so that holds 

S e (fc e ( s ))=0, 

and we assume that we can expand k e (s) using the formal asymptotic expansion 

fc e (s) = Xij + efco(s) + e 2 ki(s) + 0(e 3 ). (3.10) 

Using the expansion for S e , the periodicity of Si in y, and the Taylor series of Sq and S\ 
around (x,ko), we obtain 

S (x, fe ) + t{Si{x, k ) + k ■ W x S (x, k ) + h ■ W y S (x, k )) + 0(e 2 ) = 0. 



Tycho L. van Noorden and Adrian Muntean 



Collecting terms with the same order of e, we see that ko(s) parametrizes locally the zero 
level set of Sq: 

S (x 7 k ) = 0. 

For ki, we have the equation 

Si(x, ho) + fc • V x S (x, k Q ) + ki ■ V y S (x, k ) = 0. 
It suffices to seek for k\ that is aligned with i/q, so that we write 

fci(s) = X(s))u (s) = At 



VySo 



where, using (3.111, A is given by 

A := - 



|V„S„| 



Si _ k ■ \7 X S 

VyS \ \VyS \ 



(3.11) 



(3.12) 



(3.13) 



Each of the boundary conditions in (2.3| admits the structural form 



K(x,x/e) = for all x G T e , 



where if is a suitable linear combination of u e , Vu e , q e , p e , v e , and Vw e . Using (3.101 
and the Taylor series of K around (x, ko), we obtain 



K(x, ko) + e(k ■ V x K(x, k ) + k x ■ W y K(x, k Q )) 

+y(fcb, ki) ■ (V 2 K(x, k ))(k , h) + e 3 (...) = 0, 



(3.14) 



where T> 2 K denotes the Hessian of K w.r.t. x and y. Substituting (3.121 into (3.14), we 



can restate (3.14) in the following way: 



K(x, y) + e(y ■ V x K(x, y) + \v • V y K(x, y)) 

+y • (V 2 K(x,y))(y,\v ) + O(e 3 ) = for ally e r (a;). (3.15) 

In order to proceed further, we make use of the following technical lemmas. Their proofs 



can be found in 24 



Lemma 3.1 Let g(x,y) be a scalar function such that g(x,y) = for all y G Tq(x), 
x G f2 and t > 0. Then it holds that 

V x g = TL'^f VgSo, for x e Q, y G T (x,t). 

Lemma 3.2 Let F(x,y) be a vector valued function such that V y • F(x,y) = on 
Yq(x) := {y | Sq(x, y) > 0} and v ■ F(x, y) = on T (x) for all x € ft. Then it holds that 



r°(x) 



|V„5 | 



|V„5 | 



f • V y (^° • F)da = 0, for a; G ft. 
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Substituting the asymptotic expansions of q e and p e into (2.1 ^3), we obtain 



1 



5o 
1 



n-VyPo - nVyPx - nV x p + O(e), 



V y ■ q + V x ■ q Q + V y ■ qi + 0(e) 



0. 



(3.16) 
(3.17) 



Substituting the asymptotic expansion of q e into the boundary condition ( 2.3 5), and 
using (3.15), gives 

qo + e(qi + (V ' x qof y + \(V v q fvo) + 0{e 2 ) = 0, for all y G T (x). (3.18) 

The e _1 -term in (3.16) indicates that \7 y po = 0, so that we conclude that p is indepen- 
dent of y. Furthermore, we obtain, after collecting e°-terms from (3.16) and (3.18) and 
e -terms from (3.17), the equations for qo and p\\ 



q = -nV y p 1 - kV x Pq in Y (x), 



V y • q = 
9o = 

qo and Po y-periodic, 



in Y (x), 
on r (a;), 



where 



Y (x) :={y\S o (x,y)>0}. 



(3.19) 



(3.20) 



These equations (together with boundary conditions on the outer boundary dft) deter- 
mine the averaged velocity field given by 



q(x) 



Y (x) 



%{x,y) dy. 



Now we compute the divergence of q (where we use the e°-terms from (3.17)) 

V x Sq 



V x ■ g = V x 



Qody 



Y (x) 



■ q dy 



Y (x) 



r (x) l v a 



VyS \ 



q da 



/ V y -qidy = - vo-qi da 

Jy(x) Jr a (x) 



r (x) 
h-h, 



vo ■ ((V x qo) T y + A(V J/ g ) T ^o) da 



with 



h ■= 

h := 



vo ■ 



T {x) 



'At 



((V«<fo) 
5'i 



r (aO 



\VyS0 



£/- , , (V y g ) ^0 a CT ' 
|V y 5 | 

(Vj,qo) T fo) da. 



We apply Lemma 3.1 with g = vq ■ qo, and obtain 



Vx(vo ■ qo) = IV7 1 Vz-So, on r (a;,*). 

V w 5 
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Since q = on F (x) it follows that (V x q ) T v = ^fi^p^ Vggo, so that h = 0. 



Next we apply Lemma |3 . 2 1 with F = qo, and get consequently 

ivy o I T o • 90 - |Ty g ' Vj,(i/ • q ) da = 0. 

Again using q = on Tq(x), it follows that I 2 = 0, so that we have 

V, .g = 0. 



(3.21) 



3.4 Diffusion equation in the low diffusivity areas 



Substituting the asymptotic expansion of v e into (2.2), we obtain 



d t v = DiV y v + 0(e). 
Similarly expanding the boundary condition ( |2.3^ ), we get 

= u - Vq + O(e) on T e , 



which, after substitution into (3.15), becomes 



(3.22) 



= u — v + O(e) on r (x). 

Collecting the lowest order terms, and using that uo does not depend on y, we obtain 
the boundary condition 



vq(x, y, t) = Uq(x, t) for all y <E Tq(x), x € O. 



(3.23) 



3.5 Convection-diffusion equation in the high diffusivity area 



Substituting the asymptotic expansion of u e into (2.1 l), wc obtain 

d t u = -^D h A y u + -(Vj, ■ F h + V x ■ (D h W y uo)) 

+V y ■ (D h (V y u 2 + V x ui) -qiu - q ui) + V x ■ F h 
+0(e), 



where 



Fh ■= D h (V x u Q + VyUi) - q u . 



Using the expansions for u e , v e and v e , we first expand ( 2.3 l): 



= v c ■ (D h Vu e ) - eV • (DiVv t ) 
= -vq ■ (D h V y u ) + v ■ (D h (V x uo + V„ui)) + v\ ■ (D h V v u ) 



(3.24) 



(3.25) 



+ e(v ■ (D h (V x Ui + V y u 2 )) + v\ ■ (D h (V x u + V^iti)) + v 2 ■ (D h V y u ) - v ■ (DiV y v ) 



0(e 2 ), for all x € Pandt/ 
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Next we substitute this expansion into (3.15), and thus obtain 
= ~u ■ {D h V y u ) 

+ v ■ (D h (V x u + Vj/ui)) + ui ■ (D h VyU ) + y ■ V ' x (v ■ (D h V y u )) + Xv ■ V y (v ■ (D h V y u )) 
+ e^ ■ {D h (V x ui + VyU 2 )) + v\ ■ D h (V x u + V^tti) + v 2 ■ (D h V y u ) 

-i>o ■ (DiVyVo) + y ■ V x (vo ■ (D h (y x u a + V y ui)) + vi ■ (D h V y u )) 

+\is ■ V y (^ • (D h (y x u + VyUi)) + vi ■ (D h V y u )) 



+ ~(y,\v ) ■ (V 2 (v ■ (D h V y u )))(y,\is ) 
+ 0(e 2 ), iovyeT (x). 



Now we collect the e 2 -term from (3.24) and the e 1 -term from (3.26). Hence we obtain 
for mq the equations 



(3.26) 



A y u o = in Y (x), 

VQ-VyUo^O on r (i), 
ito w-periodic, 



(3.27) 



where Y (x) is given by (3.20). This means that u is determined up to a constant and 
does not depend on y, so that S7 v uq = 0. Collecting the e _1 terms from (3.24), the 
e°-terms from (3.26), and using that \7 y u — 0, we get for u\ the equations 



V y ■ (D h S/ y ui - q u ) = m.Y (x), 
v Q ■ (Dh(V x UQ + V^ui)) = on l (i), 
ui y-periodic. 



(3.28) 



Collecting the e°-terms from (3.24) and the e 1 -terms from (3.26), we obtain 

d t u = V y • (D h (\7 y u 2 + V x ui) - qiu - qoux) + \7 X ■ F h in Y (x), 
va ■ {D h (V x Ui + V y u 2 )) = -v x ■ (D h (W x u + Vj/Ui)) 

+v ■ (DiVyV a ) - y ■ V x (v ■ (D h (V x u + V y Ui))) 

-\v ■ V y (v ■ (D h (V x u + VyUi))) on T (x), 
u 2 y-periodic. 



(3.29) 



Integrating ( 3.29 l ) over Yq(x) and using the boundary conditions ( 3. 19 j) and (3.29>) 
yields 



\Y (x)\d t u 



V y ■ {D h (V x Ui + V v u 2 ) - qiu - q Q Ux) dy + 



Y (x) 



V x • F h dy 



Y (x) 



r (x) 



-v x ■ F h + v Q ■ (DiVyVo) - y ■ V x (vo ■ F h ) - \v ■ V y (v a ■ F h ) da 



F h dy 



Y (x) 



V^So 
T (x) \^ySo\ 



F h da. 
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Using (3.8 ), (3.13 1, and the boundary conditions (3.19 j) and (3.28 j), this can be rewritten 



as 



\Y (x)\d t u = V x • / {D h (V y ui + S7 x u ) -q u )dy 

JY (x) 

+ v - (DiVyVo) dy-h- I 2 , 



where 



h 
h 



r (x) 



y ■ V x Sq 
y ■ \7 x g - , , u ■ V y g da, 



r (x) |V„So| 



t ■ F h 



Si 



\VyS 



v -V y (i> -F h )da, 



with g := v$ ■ F^, The boundary co ndit ions ( |3 . 19| j ) and ( |3.28^ ) give us g(x, y,t) = for 
y G To(x,t). Now invoking Lemma 



integral I 2 we invoke Lemma |3.2| to 



3.1 
ob 



leads to V x g = '^^W x So. So h = 0. For the 
ain I 2 = 0. As a last step, we use the divergence 



theorem and interface condition (3.231 to obtain 



d t \Y (x)\u + / v dy = V x • / {D h (V y ui + V x u ) - q u ) dy, (3.30) 
where Y c (x) is the complement of Yq(x) in U given by Y Q c (x) := U\Yq(x) = {Sq(x) < 0}. 



Remark 3.3 Note that in this section we have not used any assumptions of the shape 
of the perforations. They may have any shape as long as their limiting shape is described 
by the level set function Sq. 



4 Upscaled equations 



The equations for lowest order terms of q e andp e , (3.19) and (3.21), v e , (3.22), u e , (3.30), 
and the coupling conditions (3.23 1 together constitute the upscaled model. In this section 
we collect these equations for the case discussed in Section[2] i.e. for circular perforations. 
For this purpose we return to a formulation in terms of r(x, t), where we use 

T (x) = {y€U\\y\=r(x)}, 
Y (x) = {yeU\\y\>r(x)}, 
Y c (x) = {y£U\\y\ <r(x)}. 



We write the solutions of equations (3.28) and (3.19) in terms of the solutions of the 
following two cell problems (see, e.g. |12|) 

A y Vj(x,y) = 
v ■ V v Vj(x,y) = -v {) ■ 
Vj(x,y) y-periodic, 



for all x G O, y G U, \y\ > r(x), 
for all iGfl, \y\ = r(x), 



(4.1) 
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and 



Wj(x,y) = V y ir j {x,y)+e j 
V v ■ Wj(x,y) = 



w 3 = 



for all x £ fi, y £ U, \y\ > r(x), 
for all x £ f2, y £ U, \y\ > r(x), 
for all x £ fi, 1 2/ ] = K 3 -)) 



(4.2) 



Wj(x, y) and irj(x,y) y-periodic 



for j = 1,2. The use of these cell problems allows us to write the results of the formal 
homogenization procedure in the form of the following distributed-microstructure model 



d t v (x, y, t) = D l A y v (x, y, t) 
d t (e(x)u + J lyl<r{x) v dy) = 



for \y\ < r(x), x £ f2, 



■ (D h A(x)V ' x u - qu ) for x £ fi, 
for i6fi, 
for x £ f2, 



(4.3) 



g = -n1C(x)W x po 
{V x -q = 

v (x,y,t) = u (x,t) for |y| = r(x), 

Uo(x,t) = Ub(x,t) for x £ T, (4-4) 
g(x, i) = qb(x, t) for x £ Y, 

J 1*0(2;; 0) = uj(x) for x G J7, 

[u (x,t/,0) = vi(x,y) for |y| < r(x), x G fi. 

where the porosity #(2;) of the medium is given by 

6{x) := l-irr 2 (x), 

while the effective diffusivity A(x) := {a>ij(x))i t j and the effective permeability JC(x) 
{kij(x))ij are defined by 



(4.5) 



and 



Oij(x) :-- 



kij (x) 



{yeU\ \y\>r(x)} 



Sij + d Vi Vj(x,y,t) dy, 



{y<£U\\y\>r(x)} 



Wji(x,y,t) dy. 



5 Analysis of upscaled equations 



In this section we investigate the solvability of the upscaled equations (4.3)-(4.5|. Note 



that the equations (4.3 ^4) for q and po, together with the boundary condition (4.4 j) are 



decoupled from the other equations. We may assume that we can solve these equations 
for q and po such that q £ i°°(S7;K 2 ) (see Assumption 2 below). Standard arguments 
form the theory of partial differential equations justify this assumption if the data q^ and 
r are suitable, see [13] for a closely related scenario. With this assumption the equations 
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( 4.3 )-( 4.5 ) reduce to the following problem 



(P){ 



9(x)dtu — \7 X ■ (D(x)X? X u — qu) 

d t v - DiA y v = 

u(x,t) = v(x,y,t) 

u(x, t) = Ub(x, t) 

u(x, 0) = ui(x) 

v(x,y,0) = v T (x,y) 



I, 



dB(x) V V 



(DiVyv) da in O, 

in B(x), 

at (x, y) € fl x dB(x), 
at i £ dSl, 
in Q, 



at (x, y) efl X B(x), 



where B{x) :— Yq(x), where Yq is defined in (3.20). Notice that in this section we again 
do not restrict ourselves to circular perforations. The perforations may have any shape 
as long as they are described by the level set 5*0. In the following sections we discuss the 
existence and uniqueness of weak solutions to problem (P). 



5.1 Functional setting and weak formulation 

For notational convenience we define the following spaces: 

Vi-H^Sl), (5.1) 

V 2 :=L 2 (n;H 2 (B(x))), (5.2) 

ffi :=£?(«), (5.3) 

H 2 :=L 2 (n;L 2 (B(x))). (5.4) 

The spaces H 2 and V 2 make sense if, for instance, we assume (like in [18]): 



Assumption 1 The function So '■ ^ x U — > R, which defines B(x) := Yq(x) in (3.20), 
and which also defines the 1-dimensional boundary Q, x dB{x) of Q, x B(x) as 

{x,y) € £1 x dB(x) if and only if S (x,y) = 0, 



is an element of C 2 (51 x U). Assume additionally that the Clarke gradient d y So(x,y) is 
regular for all choices of (x, y) £ f2 x U. 



Following the lines of [18] and 23 , Assumption 1 implies in particular that the mea- 
sures |<9.B(a:)| and |-B(x)| are bounded away from zero (uniformly in x). Consequently, 
the following direct Hilbert integrals (cf. [8] (part II, chapter 2), e.g.) 



L 2 (fl; H 1 (B(x))) := {u e L 2 (fi; L 2 {B{x))) : V y u e L 2 {n- L 2 {B{x)))} 



measurable such that 



U(X) \\h(8B( X )) < °°} 



L 2 (f2; H 1 (dB(x))) := {u : Q x dB(x) - 

are well-defined separable Hilbert spaces and, additionally, the distributed trace 
7 : L 2 (ft;H 1 (B(x))) -> L 2 {n, L 2 {dB{x))) 

given by 

~/u(x,s) := (j x U(x))(s), xeQ,sedB(x),ueL 2 (n-H 1 (B(x))) (5.5) 
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is a bounded linear operator. For each fixed x G O, the map 7^, which is arising in (5.5|), 
is the standard trace operator from H 1 (B(x) to L 2 (dB(x)). We refer the reader to 
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for more details on the construction of these spaces and to 19 for the definitions of 



their duals as well as for a less regular condition (compared to Assumption 1) allowing 
to define these spaces in the context of a certain class of anisotropic Sobolev spaces. 
Furthermore we assume 

Assumption 2 

6, D G L£(fi), 

q G L°°(Q;R d ) with V ■ q = 0, 
u b G L™(n x S)DH 1 (S;L 2 (n)), 
d t u b < a.e. (x, t) G x S, 
ui G I^(fi)nffi, 

k U Z (x, •) G L^(B(x)) n if 2 for a.e. lefl. 

We also define the following constants for later use: 

Mi := max{||u/|| L oc (n) , |K|| L °o (n) }, (5.6) 

M 2 :=max{||t;7|| L » (n ) > Mi}. (5.7) 

Note that Mi and M2 depend on the initial and boundary data, but not on the final 
time T. Let us introduce the evolution triple (V, H, V* ) , where 

V:={(£V) G Vi x 7 a I 0(a) = ^(x, y) forx £ fi, y G dB{x)}, (5.8) 
H := Hi x H 2 , (5.9) 

Denote J7 := u — u b and notice that U — at 951. 

Definition 5.1 Assume Assumptions 1 and 2. The pair (u,v), with it = U + u b and 
where (U, v) G V, is a weak solution of the problem (P) if the following identities hold 

6d t (U + u b )(j> dx + / (DV x (U + u b ) - q(U + u b )) ■ V x (j)dx = 



II v y ■ (D{V y v)(j) dadx, 

Jn JdB(x) 



dtvtp dydx 



O J B(x) 



i i DiVy ■ VyV' dydx = / / v y ■ {Di\7 y v)4> dadx 

JnJB(x) JnJdB(x) 



(5.10) 



(5.11) 



for all (</>, V>) G V and t G 5. 



As a last item in this section on the functional framework, we mention for reader's 
convenience the following lemma by Lions and Aubin [16] , which we will need later on: 



Lemma 5.2 (Lions- Aubin) Let B B Si be Banach spaces such that -Bo and B x 
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are reflexive and the embedding Bq B is compact. Fix p, q > and let 

with 

Ikllw != IN |i»(S;B ) + ll<MU'*(S;,B 1 )- 

Then W LP{S:B). 

5.2 Estimates and uniqueness 

In this section we establish the positivity and boundedness of the concentrations. Fur- 
thermore, we prove an energy inequality and ensure the uniqueness of weak solutions to 
problem (P). 

Lemma 5.3 Let Assumptions 1 and 2 be satisfied. Then any weak solution (u,v) of 
problem (P) has the following properties: 

(i) u > for a.e. igfl and for all t G S; 

(ii) v > for a.e. (x, y) £ £1 X B{x) and for all t £ S: 

(iii) u < Mi for a.e. x £ fl and for all t £ S; 

(iv) v < M 2 for a.e. (x, y) £ Q x -B(a:) and for all t £ S; 

(v) The following energy inequality holds: 

H U llL 2 (5 ; yi)nL°°(S;ffi) + \\ V \\ 2 L 2 (S;L 2 (n,V2))nL°°(S;H 2 ) 

+ \\V x u\\ 2 L 2 {s . Ml} + \\V y v\\l HSxQxB(x)) < Ci (5.12) 



where Mi and M 2 are given in (5.6) and (5.7), and where C\ is a constant independent 
of u and v. 

Proof We prove (i) and (ii) simultaneously. Similar arguments combined with corre- 
sponding suitable choices of test functions lead in a straightforward manner to (iii), (iv), 
and (v). We omit the proof details. Choosing in the weak formulation as test functions 
(<p, ip) := (— U~ , —v~) £ V, we obtain: 

\ [ ^d t U-f + \ j j d t {v-f+ j D\VU-\ 2 + { { D e \V y v-\ 2 
L Jn z JfiJB(i) Jn JnJB(x) 

= [ <i)d t u b U- + [ D\Ju b \7U- - [ V-(q(U + u b ))VU~ 
Jn Jn Jn 

< [ DVu b VU~ - j q{VU + Vu b )VU~ - [ (U + u b )dWq\7IJ- 
Jn Jn Jn 

= mmq [ |VC/-| 2 + / f7 - divgVtf- 
n Jn Jn 

- [ U + divqVU~ + / (D\7u b -u b divq)\7U-. (5.13) 
Jn Jn 
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Note that, excepting the last two terms, the right-hand side of (5.13) has the right sign. 
Assuming, additionally, a compatibility relation between the data q,Ub, for instance, 
of the type DVub — u^drvq a.e. in Q X S, makes the last term of the r.h.s. of (5.13) 
vanish. The key observation in estimating the last by one term is the fact that the sets 
{x G il : U(x) > 0} and {x G : U(x) < 0} are Lebesque measurable. This allow to 
proceed as follows: 



£/ + divqVE/ _ = 



U + divq\7U- 



{xen-.u(x)<o} 



U+divqVlJ- = 0. (5.14) 



'{xen-.u(x)>o} 

(5.15) 

After applying the inequality between the arithmetic and geometric means applied to the 
second term for the right hand-side of (5.13), the conclusion of both (i) and (ii) follows 
via the Gronwall's inequality. □ 

Proposition 5.4 (Uniqueness) Problem (P) admits at most one weak solution. 

Proof Let (it^t^), with i G {1,2}, be two distinct arbitrarily chosen weak solutions. 
Then for the pair (p, 6) := (u 2 — tti, v<i — vi) we have 



4>d t p(p 



DVpVtyS- / qpS/Lp 



dm 



n jb(x) 



D e X7 y 8V y ip = 



(5.16) 



m j b(x) 

for all (ip, ip) e V. 

Choosing now as test functions (ip, -0) ■— {p- 9) € V, we reformulate the latter identity 
as: 



-.{dtp? 



[ [ \{d t 9) 2 + I D\V P \ 2 + [ [ D e \V y 6\ 2 = [ qpVp. (5.17) 

JnJB(x) z Jn JnJB(x) Jn 



Noticing that for any e > we can find a constant c c G]0, oo[ such that 

/ q P Vp<e [ \Vp\ 2 + c t \\q\\l [ \p\": 
Jn Jn Jn 



then (5.17) yields: 



ld_ 
2~dt 



Jn Aai JnJB(x) Jn 

+ j[ D t \V y 6\ 2 <c e \\q\\l> I 
JnJB(x) Jn 



Choose 



e G 



0, min D 

nxB(x) 



(5.18) 



(5.19) 



Since for all x G Q and y G B(x) we have 9(x,y,0) — p(x,0) = 0, (5.18) together 
with (5.19) allow for the direct application of Gronwall's inequality. Consequently, the 
solutions (ui, Vi) with i G {1,2} must coincide a.e. in space and for all t G S. □ 
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Remark 5.5 At the technical level, the merit of the basic estimates enumerated in this 
section is that they are derived in the x-dependent framework and not in a fixed-domain 
formulation. Note also that the proof of uniqueness does not rely on the use of L°°- and 
positivity estimates on concentrations. 



5.3 Existence of weak solutions 

In this section, we prove existence of weak solutions of problem (P) . We will do this using 
the Schauder fixed-point argument. The operator, for which we seek a fixed point, maps 
the space L 2 (S; L 2 (0)) into itself, and consists of a composition of three other operators. 
In order to define these operators, we need the following functional framework: 

X 1 := L 2 {S;L 2 (n)), (5.20) 
X 2 •.= L i (S;H^Sl))nH 1 (S;L 3 {n)), (5.21) 
X 3 := L 2 (S; V 2 ) n H\S; L 2 (n ; L 2 (B(x)))). (5.22) 

The first operator T\ maps a / £ X% to the solution w £ X 2 of 

6d t (U + u b )(j)dx+ ( {D\7 X (U + u b ) - q(U + u b )) -\7 x <f>dx = - / f<f>dx, (5.23) 



for all (j) e fl^(n). 

The second operator T 2 maps a w £ X 2 to a solution v £ X 3 of 

/ / d t (V + w)ip dydx + / DiW y (V + w) ■ V y ijj dydx = 
Jn Jb(x) Jn Jb{x) 

u y - {DtfyiV + w))ipdadx, (5.24) 



n JdB(x) 

for all tjj £ V 2 and t £ S. 

The third operator T3 maps a v £ X 3 to / £ X\ by 



/ = / v y -V v vda. (5.25) 

JdB(x) 

The operator T : Xi X\ of which a fixed point corresponds to a weak solution op 
problem (P) is now given by 

T:=T 3 oT 2 oTi. (5.26) 
Lemma 5.6 The operator T is well-defined and continuous. 

Proof Since the auxiliary problem (obtained by fixing /) is well-posed (see e.g. chapter 
3 in |15| ), we easily see that T\ is well-defined. Furthermore, by standard arguments we 
can ensure the stability of the weak solution to the latter problem with respect to initial 
and boundary data and especially with respect to the choice of the r.h.s. /, that is T± 
maps continuously X\ into X 2 . 

Analogously, same arguments lead to the well-definedness of T 2 and to its continuity 
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from X2 to X2 C X3. The fact that the linear PDE (5.24 1 and its weak solution depend 
(continuously) on the fixed parameter x £ is not "disturbing" at this point]^] 

Since for any v £ X3 the gradient V y i> has a trace on dB(x), the well-definedness and 
continuity of T3 is ensured. □ 

Furthermore we need for the fixed-point argument that the operator T is compact. It 
is enough that one of the operators Ti, T 2 and T 3 is compact. Here we will show that T 2 
maps X2 compactly into X3. 

Lemma 5.7 (Compactness) The operator T 3 o T 2 is compact. 



Proof We will first reformulate (5.24) by mapping the x-dependent domains for the 
y-coordinate to the referential domain B(0) so that the transformed solution v is in 
L 2 (S; L 2 (n ; L 2 (B(0)))) n H l (S; L 2 (Q; L 2 (B(0)))) 

This transformation is a mapping if? : x B(0) — > il x B(x). We call 'f? a regular 
C 2 -motion if £ C 2 (Q x B(0)) with the property that for each x £ Q 

V(x, •) : B(0) -> B(x) := ^[x, B{0)) (5.27) 

is bijective, and if there exist constants c, C > such that 

c < dct\7 y ^(x,y) < C, (5.28) 

for all (x,y) £ O x B(0). The existence of such a mapping is ensured by the fact that 
S £ C 2 (tt x U), by Assumption 1. 

If 4" is a regular C 2 -motion, then the quantities 

F := V y * and J := det F (5.29) 

are continuous functions of x and y. Furthermore, we have the following calculation rules: 

VyV = F~ T VyV, 



V 

d t v = d t v, 



v y ■ j da = / JF T Vy ■ j da. 

dB{x) JT 

The transformed version of ( 5.24[ ) is now written as: let w £ X2 be given, find V £ 
L 2 {S- L 2 (n ; H^B{0)))) n H\S; L 2 (tt; L 2 (B(0)))) 

d t (V + w)ijjJ dydx+ / / JF_ 1 D i F _T V 1/ (V >r + w) ■ Vyip dydx = 
in Jb(o) Jn Jb(o) 



v y ■ {JF-*~DiF- x V y {V + w))ipdadx, (5.30) 



for all V G L 2 (n; H%(B(Q))) and t £ S. 
Denote by Tq the boundary of B(0). 



4 Note however that this ^-dependence will play a crucial role in getting (at a later stage) 
the compactness of T2. 
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Claim 5.8 F is C 2 . 



Proof of claim The conclusion of the Lemma is a straightforward consequence of the 
regularity of So, by Assumption 1. □ 

Claim 5.9 (Interior and boundary iJ 2 -regularity) Assume Assumptions 1 and 2 and take 
Vi £ L 2 (Q,H 1 {B(0))). Then 

V e L 2 (S; L 2 (Vl; Hf oc (B(0)) n Hi(B(0)))). (5.31) 
Since r is C 2 , we have 

V £ L 2 (S; L 2 (n ; H 2 (B(0)) n H^(B(0)))). (5.32) 

Proof of claim The proof idea follows closely the lines of Theorem 1 and Theorem 4 
(cf. [10], sect. 6.3) □ 



Claim 5.10 (Additional two-scale regularity) Assume the hypotheses of Lemma 5.9 
be satisfied. Then 



to 



f G L 2 (S'; J ff 1 (fi;i7 2 (S(0))nH 1 (B(0)))). (5.33) 



Proof of claim Let us take ^ CI' C Cl arbitrary such that ft :— dist(Vi' , dCt) > 0. At 
this point, we wish to show that 

V € L 2 (S; H 1 ^'; H 2 (B(Q)) n fl£(B(0)))). (5.34) 

The extension to L 2 (S; H 1 ^; H 2 (B(0)) n H$(B(0)))) can be done with help of a cutoff 
function as in flO] (see e.g. Theorem 1 in sect. 6.3). We omit this step here and refer the 
reader to loc. cit. for more details on the way the cutoff enters the estimates. To simplify 
the writing of this proof, instead of V (and other functions derived from V) we write 
V (without the hat). Furthermore, since here we focus on the regularity w.r.t. x of the 
involved functions, we omit to indicate the dependence of U on t and of V on y and t. For 
alH G S, x G O' and Y e Yq, we denote by and Vfi the following difference quotients 
with respect to the variable x: 

U h (X,l) .— 



ft 

V"(a; + he i: y,t) - V(a;,j/,i) 



ft 

We have for all ^ € £ 2 (ri', #o (-B(O))) the following identities: 

J(x + hei)d t (V{x + ha) + U(x + fte,-))^ + / 5(x + fte l )V y y(a; + he^Vyip 

Q'xB(O) JSl'xB(O) 

/ ^ • (S , (x + fte 4 )L> £ V y V r (a; + fte l ))i/;d ( T = (5.35) 
'n'xr 
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and 

/ J(x)d t (V(x) + U(x))^+ [ S{x)V y V{x)V y i> 

Jn>xB(0) JffxB(O) 

- / v y - (S(x)D e \7yV(x))tpdcr = 0. (5.36) 

Subtracting the latter two equations, dividing the result by h > and choosing then as 
test function ip := V£ yields the expression 

A 1 +A 2 + A 3 = 0, 

where 

A 1 := f VI [J(x + ha)dt(V{x + ha) + U(x + he,)) - J(x)d t (V(x) + U(x))] \ 

JO'xB(O) ft 

= / Vi(d t Vi + d t Ul)J(x)+ [ (d t V(x + h ei ) + d t U(x + hei))4(x)Vi 

JO'xS(O) JffxB(O) 

A 2 := [ hs(x + hei)VyV(x + kei)-S(x)VyV(x)]VyVi 

Jn>xB(o) h 



Jn'xB(o) Jn'xB(o) 

A 3 :=- [ \v y ■ [S(x + he t )V y V(x + he t ) - S(x)V v V(x)] V l h 
Jn>xr h 

= I vy [Sy7 y V{x + he,) + SVyViVl). 

Re-arranging conveniently the terms, we obtain the following inequality: 

\ I iV?f\J{x)\ + / \S{x)\{V y Vff < [ \v£d t uiJ(x)\ 
L Jn'xB(o) Jn>xB(o) Jn>xB(o) 

+ [ \{d t V(x + he i ) + d t U(x + he i ))Ji{x)V^\ 

Jn'xB(o) 

+ l \SiVyV(X + h ei )VyVi\ 

JS)'xB(O) 

+ / \vy (SVyVi)Vi\ + I \vy {Sy7 y V{x + hei)VZ)\ 

5 

= J>. (5-37) 

£=1 

To estimate the terms If we make use of Cauchy-Schwarz and Young inequalities, the 
inequality between the arithmetic and geometric means, and of the trace inequality. We 
get 

|2 



\h\ < ( " xB(0)) ||Vfe|l L 2(n/ X B(o)) + ^\\dtU l h \\ L ^{n' xb{o)), (5-38) 
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II J||2 

\h\ < £ °°^' x5 (°» 2 {\\atV(x + kei)\\^ in , xBm + \\d t U(x + he i )\\ L 2 {n , xBm ) 

+ II^IU»(n'xfl(o)), (5-39) 
|/s| < e||V y V^|||2 (n , XjB(0)) + c e \\Sl\\l x{Q , xBm \\V y V(x + hei)\\ 2 L2(Q , xBm , (5.40) 

Wy ■ (SV y Vi)V£\ < ||5|| L - ( n'xro)II^IU-(n'xr ) / Wv • V„^| 

Q'xr JSl'xr 

< |B(0)|5||S'|| i oo( n , xro) ||V h l || i oo( f2 , xro) ||V h l ||| 1(n , ;iJ - 2(B(0))) , 

(5.41) 

and 

= (f2'xro)ll^!llL=(a'xr )ll^llii(0';H 2 (B(o)))- 

SI'xTq 

(5.42) 

Note that all terms are bounded from above. To get their boundedness we essentially 
rely on the energy estimates for V, U, as well as on the L°°-estimates on V and 
V/l on sets like SI' x 5(0) and U' x T . The conclusion of this proof follows by applying 
Gronwall's inequality. □ 



Using the claims above, we are now able to finish the proof of Lemma |5.7[ by noting 
that T 3 oT 2 : L 2 {S;H l {il;H 2 n Hq(Bq))) -> L 2 (S;ff x (n)) is continuous and compact 
via applying Lemma 5.2 with Bq — H 1 ^) and B = B\ = L 2 (il). □ 

Putting now together the above results, we are able to formulate the main result of 
section [5j namely: 

Theorem 5.11 Problem (P) admits at least a global- in-time weak solution in the sense 
of Definition \5.1\ 



6 Discussion 

The remaining challenge is to make the asymptotic homogenization step (the passage 
e — > 0) rigorous. Due to the x-dependence of the microstructure the existing rigorous 
ways of passing to the limit seem to fail [3j [14 21 . As next step, we hope to be able 
to marry succesfully the philosophy of the corrector estimates analysis by Chechkin and 
Piatnitski [6] with the intimate two-scale structure of our model. 
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